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Abstract. The Dirichlet problem with oscillating boundary data 
is the subject of study in this paper. It turns out that due to integral 
representation of such problems we can reduce the study to the 
case of surface integrals of rapidly oscillating functions, and their 
limit behavior: 




where g(x, y), represents the boundary value in the Dirichlet prob- 
lem. The lower dimensional character of the surface T produces 
unexpected and surprising effective limits. In general, the limit 
of the integral depends strongly on the sequence e = £j chosen. 
Notwithstanding this, when the surface does not have flat por- 
tions with rational directions a full averaging takes place and we 
obtain a unique effective limit in the above integral. 

The results here are connected to recent works of D. Gerard- 
Varet and N. Masmoudi, where they study this problem in combi- 
nation with homogenization of the operator, in convex domains. 
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1. Introduction 

1.1. Background. In this paper we consider Dirichlet and related 
problems, with rapidly oscillating data. Although we treat the case of 
Laplace equation, our analysis straightforwardly extends to general 
cases for equations that admit a Poisson/Green representation. Such 
integral representation, in turn, reduces the study of the problem to 
the corresponding integral equation, and hence the analysis of the 
integral of rapidly oscillating functions becomes central. 

To fix ideas, let T be a C 1 surface in W 1 (n > 2), not necessarily 
bounded. We assume g(x, y) is integrable in both variables over f, 
and 1-periodic in y-variable i.e., g{x, y + k) = g(x, y) for k e Z n . In this 
paper, we shall study possible limit behaviors of the integral 

a) 5s X* ( y/ f) d(Ty/ 

and we shall prove that under mild conditions on the surface T there 
is an effective limit as e tends to zero. 

In general we show that the above integrals stay within the bounds 

of the interval [ j T gXV' v y)> j T g*(V' v y)\ where ~g are defined below, 
as the infimum respectively supremum of the average-integrals of 
g(y, -) over closed loops of the plane {x : x • v y - 0} on the torus; here 
v y denotes the normal vector of T at y. 

To illuminate the application of this to the Dirichlet problem, let 
us consider a bounded C 1 -domain Q c W (n > 2). Let further g(x, y) 
and f(x, y) be continuous in both variables, and periodic in y. Let u e 
be a solution to the Dirichlet problem 

(P £ ) Au £ = -\i £ in Q, u e = g(x,x/E) on <9Q, 

where \i £ - f(x / x/e)xT do x/ and T is a C 1 surface compactly in Q. 
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Since solutions to this problem can be represented by surface in- 
tegrals of g and / (through Poisson and Green functions), we may 
directly apply our results. Hence, a by-product of our surface inte- 
gral homogenization, is nontrivial limit scenarios, as e tends to zero, 
of equation ( [P^] ) and its solution. 

It is noteworthy that, through integral equations of Fredholm type, 
or just standard functional minimization, our technique applies to 
homogenizations with oscillating Neumann data (see Section |63| ). 
For fully nonlinear operator the Neumann problem has been treated 
in |CKL| . see also IIBDLSp . 

Another issue that occurs along such analysis is the study of the 
speed of convergence. In problems, where governing partial differ- 
ential equations have oscillating coefficients, one uses the standard 
method of expansion 



in an appropriate space. Here Uu, e is the so-called boundary layer 
term, and finding it is part of the problem. The function Uu, e will then 
solve a Dirichlet problem with oscillating boundary data. We refer 
to (AA) for some backgrounds and to [GMJ for recent developments 
in the topic. 

Remark 1.1. It is noteworthy that our approach applies to general equations 
of the type 



with oscillating Dirichlet/Neumann data. This, even though straightfor- 
ward, becomes quite technical and is therefore outside the scope of this 
paper. Thus for clarity of the exposition we shall treat the Laplacian case, 
only. 

1.2. Heuristics. Let us set x = y(mod 1) if x e = {0 < X\ < 1}, and 
y - x € Z n , and adopt the standard notation from homogenization, 
as well as that of ergodic theory. 

For clarity of the ideas, let us also deal with problem ( |P^| ) in the 
case / = 0, and see how it reduces to the integral averaging, and then 
how the limit integral is obtained. 



Since harmonic functions can be represented using Poisson kernel 
P(x,y), 



we must then analyze the behavior of this integrals as e tends to zero. 



u e (x) = Uq + eu\{x,xfe) + euu, e + 0(e), 



div(A(x/e)Vu) = f(x/e) 
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For x £ <9Q, P(x, y) is continuous and we can rewrite the integral 

(2) ue(x) » Y P(x, yJ)\dQ n Q (yOl f g(y>, y/e) Ar y/ 

y JdClnQrftj) 

where r ; is small enough, independent of a, and Q r .(y ; ) is a cube with 
size r 7 and center y ; . This obviously brings us to integral ([!]). Observe 
that at this stage we cannot replace the second variable of g with yi, 
due to rapid oscillation for e « 0. 

Further assume that <9Q n Q r .(y ; ") is so flat that at e-scale it is ap- 
proximately e 2 -away from its tangent plane; this requires a C 2 -graph 
locally, but our techniques/proofs work for C 1 -surfaces. 

The idea now is to cover the boundary dQ n Q r .(y ; ) with finite 
many, and small enough cubes, so that each part of the surfaces is as 
flat as we want to. In particular we have 

f g(y j , yle) doy = Yf g{y j , y/e) do y , 

JdQnQ rj {yi) i JdQnQ e . r . (y^) 

where E\ is to be chosen small enough. This part is slightly more 
delicate, and needs extra care. We would also prefer to take Ei - V^- 
The trouble is that this might not work, as the points are changing, 
and we may loose lots of information about the behavior of the 
surface dQ n Q £ . ry (y z,; ). Indeed, an scaling of the integral gives 

JdnnQ (e . /e)r .(y l >0 

and the question is whether this integral will converge to the mean 
f + g(yi, y)dy. This would happen exactly when the normal vector of 

the surface dQ at y z,; is irrational. In other words the surface <9Q n 
Q(£i/E)rj(y l '0 (which is almost a plane) will foliate the n-dimensional 
torus, as e tends to zero, provided we have chosen Si/e — > oo, and the 
normal vector of the surface <9Q at y z,; is irrational. In particular 

(3) lim + g{yKy)do y = ^(y ; ). 

^ u J^QnQ (£ . /£)r .(yW) 

If the set of points on dQ with rational direction has zero surface 
measure (this may fail if there are so-called flat spots with rational 
directions) then we could cover the boundary with small cubes cen- 
tered at points with irrational normals. From here combining 
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we obtain 



It is apparent that if there is a flat portion of the boundary with ra- 
tional normal, then a full foliation cannot take place at such portions. 
Consequently, the resulting (mod 1) surface will be a close simple 
curve over the n-dimensional torus. Hence different sequence of e 
may give different shifts of this close loop, and hence the possibility 
of a parameter family of values, in between [^(y, v y ), ~g(y, v y )]. 

Remark 1.2. A word of caution: As e tends to zero, one may obviously 
rescale the integral, by a change of variables, as we did in our heuristic 
explanations above. This scaling makes the surface to be scaled (we assume 
the origin is NOT on the surface) so as to disappear in the limit. This 
naturally would make it impossible to compute the limit integral. However, 
the periodicity of the function g in its second variable, implies that we can 
bring back the surface so as it passes through the fundamental cube using 
(modi) argument. In particular this means that despite the integration is 
on a fixed surface, the surface itself will start jumping forth and back due to 
the variable e. 

It would be a good idea for the reader to consider simple examples such as 
integration over a line segment in the plane, by varying the normal direction 
of the plan from rational to irrational. 

1.3. Plan of the paper. In the next section we shall introduce all def- 
initions and notations. We take care of some technicalities in Section 
[3j We shall formulate our main results concerning surface integration 
in Section |4j and its generalization to other type of functions, such as 
layered-densities, almost periodic functions will appear in Section 5 
Several interesting applications to PDE are mentioned in Section 6 
and several Examples are also given in Section [7| 



2.1. Averaging and Ergodic theory. 

Definition 2.1. Let vbea vector in W, z a fixed given point, and Q r (z) = 
{x : \Xj - Zj\ < r}. Let g(x, y), be integrable in both variables over the plain 
{(y - z) • v = 0}, and periodic in y-variable. We define 




2. Definitions and Notation 




and 
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Later we shall consider cases where v - v z , is the normal vector at z on a 
given surface T. We also define the average of g as 

gte) = f g( z >y) d v> 

forQ\ = {x : < x { < 1}. 

Henceforth we shall assume all vectors have length one, unless 
otherwise stated. We shall also without loss of generality assume 
that the surface is orientable, and fix a consistent choice of normal in 
clockwise direction (this choice is obvious if T is the boundary of a 
domain). 

It should be noted that if the set {(y - z) • v - 0}(mod 1) foliates the 
cube then the integral 

lim \ g(z,-\do y 

r ^°° J{(y- Z yv=0}nQ r (z) V £ ' 

converges to the average g(z), and this happens exactly when v is 
irrational direction (see Lemma |3.2| ). When the set {(y - z) • v = 
0}(mod 1) does not foliate the cube Q| (or the n-dimensional torus) 
then we shall get the limit as an integral over a closed loop flow on 
the torus. In particular the value of the integral exists, and depends 
on v and e. For different values of e this loop translates over the torus 
and will give rise to supremum respectively infimum values of the 
integrals, as defined above by ~g , g^ respectively. 

Remark 2.1. It should be remarked that in the definitions of the above 
averages ~g , g^ we could have replaced 



f 



{(y-z)-v=0}nQ r (z) 



by 



f g( z >^) d °y> 

J {(y-z)-v=0}nD v L 7 

for any domain D containing the point z. The reason for this is that the set 
T = {(y - z) • v = Q)(mod 1) either foliates the whole n-dimensional torus, 
or it is a closed simple curve on the torus. In either case, when e tends to 
zero, the piece of plane {(y - z) • v = 0} n D will have the same effect as T. 

Since the direction of the normal of the plane {(y - z) • v = 0} will 
play a crucial role in our analysis, we introduce proper definitions 
well-known in Ergodic theory. For readers 7 convenience we also 
prove some of these well-known results here. 
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Definition 2.2. (IDDC) 

(1) The vector v e S n_1 is a rational direction if there is m e Z n such 
thatv=^. 

(2) The vector v e S n 1 is an irrational direction ifv is not a rational 
direction. 

(3) A smooth surface T satisfies Irrational Direction Dense Condition 
(IDDC) if 

o T (T\{x e T :v x is a irrational direction. }) = 
for the surface measure o T . Here v x denotes the normal to T at x. 



3. Technical preliminaries 

In this section we shall recall some standard facts from Ergodic 
theory, and also state and prove some averaging results that will be 
needed for the proof of the main theorems. 

The first lemma is a version of WeyTs Lemma and is well-known 
fact about uniform distribution. 

Lemma 3.1 (WeyTs Lemma). Let v = (vi, • • • , v&) e (d > \), he such 
that (l/Vi,-- • ,Vd) is an irrational direction. Then, for h e L 2 (0,1), and 
a k - Zf=i kiVi(mod 1) we have 

— V h(a k ) I h(t)dt, k = (k lf ■■■ ,k d ) 
M ££. Jo 



asN ^ co. HereI N = [k e Z d : |Jt/| < N,z = 1,- • • 

That (1, Vi, • • • , Vd) is an irrational direction, means that at least one 
of the numbers v ; is irrational. In this paper d = n - 1 is the only case 
that is used. 



Proof. Set 



to be the sum of Dirac delta functions. First, let us consider the case 
h — h m - e 2mmt and m ± 0. By the assumption there is v\, which is an 
irrational number, and hence 



(4) |MW|^nJ =1 



N 



kj=-N 



inmkjVj 



< 



c 



2N\1 - e 2inmv i\ 
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since mv\ cannot be an integer i.e. 1 - e 2mmvi ± 0. Now by Fourier 
expansion (extending h as 1-periodic function), for some b m e IR, 

h(s)= f h(t)dt + Yb m e 2inms , 

which in combination with Q results in 

h(t)dt + Yb m ^ V/N ((? 2nms )^ h(t)dt. 

m^l J 



Lemma 3.2. Lrf g, g , g^ and g be as in Definition 2.1 Then the following 
hold: 

(i) Ifv x is an irrational direction, then 

f(x/v x ) = gXx,v x ) = g(x). 

(ii) The following inequalities always hold 

~g (x,v x ) >g(x) >gX*,v x ). 
Proof. First note that we can replace the integral in the definition of 

&/ and f b y 

g^X^dOy, 



j 

J\(v 



!{(y-xyv=0}nQ N 

with N positive integers, and N — > oo. This simplifies the matter 
slightly. 

Set n = n(v, x) - {(y - x) • v = 0}. We may without loss of generality 
assume v = (1, V). Since v is irrational direction, Tl(mod 1) will foliate 
the unit cell, according to standard foliation theory and results in 
integral average over Q+ . To see this we assume (by using periodicity) 
that x e Q+ and that the plane IT cuts the Xi-axis (or some other axis), 
let a be the point of intersection between the Xi-axis and this plane, 
so that n(v,x) - {x : X\-a^-V- x'\, and (l,v') is irrational (by the 
assumption), and < a < 1 due to periodicity of g(x, y) in y. Having 
fixed v, we shall now use the notation Tl t = IT(v, te\), the plane with 
normal v through the point (£, 0'). 

For k! e define Sjf to be the cylinder generated by the Xi-axis 
and the (n - l)-dimensional cube k! + {y : < y,- < 1, z = 2, • • • , n} 
(these are cubes in translated by the vector k'). Set further 
a k , - a - k f - v'(mod 1) (k f = (fc 2 , ■ * * ,^n))/ then n % , = (Ft n S^) (mod 1). 
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Let N > be a large number, and I' N = {k' e Z n_1 : \k\ < N}. Then 
by periodicity of g(x, •) 



I g(x,y)do y = V I g(x,y)da 

Set now K;(a^) = w(x, a^) := L g (x, y) da y . Then 

JV W> 

f* 11 % ^ 

jQ N nn \ v \ UatI p^T 



N 



where we have used \Q N n n| = (2N) n_1 |v| = I^IM. 

Since v is irrational, from Lemma [3TT1 we conclude 
(5) 

lim \ g(x,y)do y = — I w(t)dt = — I I g(x, y)da y d£ = g(x), 
N ^°° JQ N nn * l v l Jo M Jo Jn, 

where 1% = I^nSc (modi). Therefore lim^ f [iy . x) . v=0]nQr g { x > V lj da y = 
g(x) which is independent of e. Now we have conclusion (i) from the 
definition of ~g(x, v) and ~gX x > v )- 

To prove (ii) let us suppose that v is a rational direction, otherwise 
the conclusion follows by the equality in (i). Since v is rational, it is 
not hard to see that the restriction of g(x, y) on this hyperplane will 
be periodic with period T, say; for g(x,y/e) the period will be eT. 
In particular the integral can be seen as integration over a spiral-like 
plane on the torus (a closed loop). The two dimensional case can be 
illustrated by a curve on the torus that loops over itself. In particular 
the limit integral (w.r.t. r) in the definition of g^ ~g can be replaced by 

gX x > v x) = lim T g U, ~ ) do y = lim jr g (x, y) da y/ 

e J{(y-x).v=0}nQ £T \ Z' ^ J{(y-x £ )-v=0}nQ T 

g (x, v x ) = lim -T g (x, -)d(jy = lim -f g (x, y) da y . 

£ J{(y-x)-v=0}nQ £T V £/ £ J{(y-x £ )-v=0}nQ r 

An observation here is that the supremum value of the mean integral 
w.r.t. e actually is taken for some value er 0/ and naturally for many 
other values, due to periodicity in e. One can think of the situation 
as parallel planes in IR n with fixed distance e\ from each other are 
moving, simultaneously by keeping a fixed distance between them, 
in the orthogonal direction of the plane (mod 1), when e changes. 
Obviously, these planes foliate once e ranges [0, 1). 
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From here one deduces 



Jq+ Jo Jn(v,x e )(mod i) 

< f g(x,y)do y <g\x,v). 

Jll{v,x eQ ){mod 1) 

Similarly 

f g{*ry)dy>gX x rV). 



In our analysis of the limit behavior of the integral ([!]) we will use 
Definition |2.1 1 in a slightly different way Next lemma will give us a 
hint in that direction. 

Lemma 3.3. Let g be as before, n(v, z) = {(y - z) • v = 0}, and z eT. Then 
for any R £ /* oo (as e \ 0) we have 



limsup \ 



g(z,y/£)do y < g{z,v), 

.z)nQ £Re (z) 



and 



lim inf \ g{z, y/e)da y > g # (z, v). 

Jn(v,z)nQ £Re (z) 



In case v is irrational we obtain 



lim-f g(z,y/e)do y = g(z). 

Jn(v,z)nQ ERE (z) 

We leave the the reader to verify this simple fact. 

We shall next make the previous lemma even more general by 
letting the plane be replaced by very smooth surfaces. Let T be a 
smooth surface, with module of continuity t = Tr for its C 1 norm. 
Define 



/ 1 1 

(6) M e = A /min(— -, -) / oo, as e \ 0. 
Then 

(7) M e z(eM e ) < ^t(V^) 0, as e 0. 

Lemma 3.4. Let g be as before, and T a smooth C 1 surface, with module of 
continuity t = x r /or its C 1 norm. Let further p £ = eM € where M c is as in 
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([6]). Then, for z e T, and t] > there exists e Vz>n such that for all e < e Vz/7] 
we have 

gX z > v z )-tj< -If g (z, -)do y < g\z, v z ) + 77. 

JTnQ pe (z) V Z> 

Proof Set z e - - £ (mod 1), and T e z :- {x : e(x + z) e T}. Then we have 

7 g\ z >-) da y= T g ( z > z ? + y) d(J y/ 

JrnQ Pe (z) V JnnQMeCO) 

where 

r 2 £ n Qm £ (0) c jx € Q Me (0) : |z • v 2 | < = M,t(^M,)| . 

If we set n Vz = {x : x • v z = 0}, then by continuity of g(z, .), and that for 
e small enough the surface T £ z is M £ t(^M £ )(^ 0) close to Tl Vz in B Me (0) 
we will have 

r g(z,z £ + y)da y = -T #(z,z £ + y)da y + 7, 

JrfnQ Me (0) ^n V2 nQ Me (0) 

where 

1 = f (S"(z,z £ + y)-g(z / z £ + y))da y . 

JrfnQ Me (0) 

Here y = y + v z |y|s £ with s £ = t(^M £ ). Estimating I we obtain 
|I| < Cz g (\y\s £ ) < Cz g (M £ s £ ) < 7] 1 2, for e small enough. Here z g is 
the module of continuity for g From here, and Lemma [33] the state- 
ments in the lemma follows, provided we have taken e small enough 
depending on v z . □ 

4. Surface integrals of oscillating functions 

Theorem 4.1. Let T be a C 1 surface in IR n , and g(x, y) be integrable in 
x-variable over T, and continuos and 1-periodic in y in R n . T/zen 

limsup | g[y,-\doy< I g*(y,v y )da y , 

and 

Moreover if! satisfies IDDC then an effective limit exists and we have 
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Proof. We shall prove the limit superior estimate only. The limit 
inferior estimates follows in a similar way. The last statement follows 
in an obvious manner. 

Let us fix a small positive constant r\ > 0, to be decided later. 
Without loss of generality, we may assume T is bounded, and g(x, y) 
is uniformly continuous function on T X IR n . 

Next from Lemma |3.4| for each z e T we may take £ Vz/7] such that 
the estimate 

(8) gXz, v z ) - T] < \ g [z, 7 W y < g\z, V z ) + 7]. 

JrnQ Pe (z) v 

holds for all e < e Vz/7] . Next we cover T by cubes Q Pe (z) for all z e T 
and all p e - eM £ with e < e Vzr7] . This is a fine cover of T, and therefore 
there is a finite disjoint sub-cover such that T \ U^ x Qp (z ; ) has surface 

measure less than 77. Here z ; e T, and 8j - e v .. Next 
and by ^ we obtain 

N 

g(y,y/^)da y < ^ (gV,v z ;) + r/) |T n Q Pe .(z j )\ + r/HglU 



< 



J" ^"(y/ v y)^y + (1 + 1 1^1 1- + |r|)r] < J" g\y,v y )do y + rj , 



for any rj as small as we wish, provided 77 is small enough. 
In a similar way we obtain the estimate from below 

As r/o was arbitrary we have the two main estimates in the statement 
of the theorem. 

Putting these together along with the IDDC we shall have the 
third statement. Indeed, due to IDDC we have that the set jieT: 
v x rational } has zero surface measure, and the integral over this set 
is zero. For the rest of T we have irrational normals only, and hence 
the full averaging (Lemma 3.2 ) takes place and we obtain g^ - g*. □ 
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5. The case of Layered densities, almost periodicity, and 

ergodicity 

In this section we shall deduce, similar results of that in The- 
orem |4.1| while replacing the periodicity assumption with layered 
materials/densities, almost-periodic and ergodic case. Nevertheless 
we shall only mention the results without deepening much into the 
analysis. The reader may easily verify the statements. 

5.1. Layered Materials. If we assume the function g{x,y) is inde- 
pendent of (i/fc+i, • • • , y n ) and is 1-periodic in (yi, • • • ,y k ) then one 
may naturally obtain results reminiscent of that of layered materials 
in homogenizations for PDE. Indeed one can obtain the following 
obvious result: If the surface T does not have any flat parts in directions 
e { (i = 1, • • • , k), then the averaging takes place. 

5.2. Almost Periodic Case. In the case of almost periodic functions 
(see [Bohr] ) one obtains similar results by replacing the average inte- 



gral g(x) with 

g(y) = lim + g(x / y)dy. 

The obvious details are left to the reader. 

5.3. Ergodic Case. The results of our main theorem can be gener- 
alized further to the case of functions with ergodic properties (see 
BJKOfl ). Indeed, if we assume g(x, y, co) is defined on T x W 1 x D, for 
some D cW and that g(x, y, co) is statistically homogeneous field in 
(y, o;)-variable, i.e., g(x, y, co) - h(x, T y co) for some random variable 
h(x, z) (random w.r.t. second variable) with underlying probability 
space and an ergodic n-dimensional dynamical system T y . Hence 
h{x, Tyco) admits an averaging 

lim *r h(x, T y / £ co)dx - Exp(/z(x, •))• 



One may now deduce similar results as that in Theorem 4.1 with 

g(z) = Exp(£(z, -,.)). 



6. Applications to Partial Differential Equations 

6.1. Dirichlet problem: Elliptic Case. Let Q be a bounded domain 
in R n , with piecewise smooth boundary. Let g{x, y) be as before, and 
f{x, y) have the same property as g. For simplicity, we shall assume 
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/, g are continuous in both variables (actually L 2 (<9Q) would suffice). 
Let further r c Q be a piece wise C 1 -curve, and define 

V* =7 Xr do, /J„ = J \XyA° J 1 = lxi,d°' 
Then the following result hold. 
Theorem 6.1. For a solution u £ of §P^ , we define (in Q) 

u*(x) := limsupM e (x) and u*(x) := liminf w e (x) . 

£ ->o 

T/zen the following hold (in the weak sense): 

Aif > -jS", At/" < in Q, 

u*(x) < ~g(x,v x ) on <9Q, u*(x) > gX x ' v x) on <2Q. 

Moreover for any sequence of u e , there is a subsequence converging to a 
function u in Q, satisfying 

-Jl < Au < -/I in the weak sense in Q, 
g*( x ' v x) ^ lim m(z) < limu(z) < ~g(x,v x ), x e dd. 
Proof By Green's and Poisson's representations we have 

Ue(x) = p ( x ' y)g (]// \ ) do y + G(x, y)f (y, V - j da y , 

where P and G are Poisson respectively the Green functions of the 
domain. 



By Theorem 4.1 



u\x)< I P(x,y)g\y,Vy)do y + \ G(x,y)f (y / v y )do y , 
Jdn Jr 

and 

u*(x)> I P(x,y)gXy,Vy)do y + I G(x,y)JXy,Vy)do y , 

which implies the first two statements in the theorem. 

Now the last statement follows from the above inequalities, in an 
obvious way. □ 

To obtain an effective limit, in the above theorem, one needs to 
consider domains satisfying IDDC. In particular, using the above 
theorem and Lemma [372] (ii) we have the following result. 
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Corollary 6.2. Let dQ and T satisfy IDDC, and u e be a solution to problem 
( |P^| ). Then, u E converges to a function u in Q, satisfying 

Au = -Ji(x), in Q, u = on <9Q. 

The general nature of the method employed here suggests that 
we can apply this to situations where integral representations are 
possible. This can naturally go beyond the Dirichlet problem, or 
the Laplace operator, and as general as to systems, and equations 
of higher orders. We state explicitly that once one has an integral 
representation of any function, through kernel functions, then one 



may conclude similar statement as that of Theorem 6.1 and Corollary 



6.2 We leave it to the reader to apply this to their favorite scenarios. 



6.2. Dirichlet problem: Parabolic Case. Let us now consider the 
case of parabolic equations 

(9) Au E - d t u E - -\i E (x, t) in Q, u{x, t) - g E (x, t), 

where d\i E {x,t) = f E {x,t)xi Q do, f E {x f t) = f(x,t,x/e,t/e), g e (x,t) = 
g(x, t,x/e,t/ e), and f{x, y, s, f), g(x, y, s, t) are integrable in {x, s)-variables 
over T and dQ respectively, and 1-periodic in (y, £)-variable. 

Here one may consider different cases, such as Q = D x (0, T), with 
D a domain in R n , or Q c R" x R and time varying. Also T can be 
take to be either time independent or varying in time. 

Now a similar argument, using Poisson and Green representa- 
tions, as in Section [6J] can give us various type of results. In the case 
Q = D x (0, T) one obtains same type of results as that of Theorem 
6.1 To obtain results of the nature of Corollary |6.2| one needs 



• either to assume Q = D x (0, T), and D has IDDC condition 
along with /, and g being independent of their fourth variable, 
i.e. f e (x,t) = f(x,x/e,t), g E (x,t) = g(x,x/e,t) 

• or Q and r , both have IDDC condition in 

In addition, Corollary |6.2| can be obtained similarly for f e (x,t) - 
f{x, x/e, t/e 2 ), g E (x, t) = g(x, x/e, t/e 2 ) with Q = D x (0, T), and D has 
IDDC condition. 



6.3. Neumann problem. As mentioned at the end of the previous 
section, one can apply the technique of this paper to far reaching sce- 
narios and problems, involving integral representations. The Neu- 
mann problem naturally fits into this category through Fredholm's 
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alternative, and an integral representation. Indeed, let u £ be a solu- 
tion to the problem 

Au £ - 0, in Q, d v u £ - g(x,x/e) on <9Q, 

with g satisfying compatibility condition j dQ g(y / y/£)da y - for all 
e > 0. Then it is well-known that 



u e (x)= | F{x,y)(p £ {y)do y , 
JdQ 



where F is the Fundamental solution for Laplace operator (or the 
corresponding operator), and <ft solves the Voltera integral equation 
of second kind, i.e., 

(f> e (x) = I d v F(x, y)(p £ {y)da y + g(x,x/e). 

JdQ 

As e tends to zero, <p £ tends (weakly in L^^Q)) to a limit <p solving 

(j) (x) = I d v F(x, y)(po(y)do y + g(x), 

in a weak sense over the boundary of Q. This happens exactly when 
the boundary of Q has IDDC. 

More accurately, the kernel of the bounded operator 

T{v) = v(x) - I d v F(x, y)v{y)do y - g(x, x/e), 

JdQ 

acting on I}(dQ) space, is upper semi-continuous and has a unique 
element. In particular lim £ ker(T £ ) c ker(T°), where 

T°(v) = v(x) - I d v F(x, y)v{y)do y - g{x). 

JdQ 

By uniqueness of the solutions to the Fredholm operator this kernel 
must have only one element, and hence <ft £ — > <po, with <ft e Ker(T°). 
In other words <fto solves the Voltera equation above for the function 

From here it follows that u £ converges to Uq = f dQ F(x,y)<fto(y)d0y, 
with <po solving 

cj) (x) = I d v F(x, y)(j) (y)do y + g(x). 

JdQ 

Hence Uo solves the averaged/effective Neumann problem 
Auq = 0, in Q, d v u £ = g(x) on <9Q. 
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A different way of analyzing this is to consider the solution of the 
Neumann problem in the weak form 

I Vu £ (y)-V(p(y)dy= | g(y / y/e)(p(y)da y/ 

Jo J<9Q 
where <p is a test function in a reasonable class. Letting e — > we see 
that the integrals converge to 

I Vu (y)-V(p(y)dy= | g(y)(p(y)do r 

Jq J<9Q 
The latter in turn solves the Neuman problem with g(y) as the amount 
of flux at each boundary point. 

7. Examples and illustrations 



The behavior of the limit integrals in 2.1 are directly related to 
foliation of the fundamental cell Q+. To illustrate this (in R 2 ) consider 

a sequence pj - (- V2//, 1). For a e [0, 1), let T. be the line through 
the point (0,0) and orthogonal to pj. Then, due to the fact that pj is 
rationally independent Lemma [372] implies 

lim \ g(x,y/e)do y = g(x). 

JlpQ P£ 

On the other hand Z ; — > Z which is a line through (0, a) and parallel 
to the Xi-axis. For the limit of the average for this line we then have 

lim \ g{x, y/e)do y = g(x, a). 

^° Jl nQ P£ 

where averaging takes place only on y x -variable. In general ~g{x, a) ^ 
g(x). More importantly, we can in general not expect to take a very 
small £ and expect 

f g(x,y/£ )do y <g\x), 

for estimates from above, or > gX x )' when considering estimates 
from below. The reader may easily verify this by taking g(x, y) - 
g(y) = | sin nyi sin ny 2 \, and make computations to arrive at 

lim -r | sininyx/e) sm(ny 2 /e)\da y = 

and 

C . . / / x . / , mj 2sin(raz) 

lim -f | sm(7Tyi/6')sin(7iy2/6 , )|a(J v = . 

-'/onQp, n 
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Let us set 

L a j = I) n Q Peo/eo (mod 1). 
Then one readily verifies that for £ fixed, and very small 

L a .cQ+n\x:a <x 2 <a + -A 1 

1 \ J 1 

where M £Q - p £o /£o- In particular, for j large L ; will never foliate the 
unit cell, and hence it is impossible to approximate the integral over 
SO by any covering, however small. 

7.1. Example 1. We consider the case when g is periodic only in X\- 
direction and when the domain is a slab with a unit normal direction 
v. 

For ave S n_1 , set Q = {x : -R 1 < x • v < R 2 }. Let g(x) = g{x x ) be 
independent of (x 2 , • • • ,x n ) and 1-periodic, i.e. g(x + k) = g(x\ + k\) - 
g(%i) = g(x) for hZ n . Now let u £ be a solution of the following 
equation 



(10) 



Au £ - in Q, 

u, (x) = M on x • v = -R {/ 
u e(x) = g(f) onx- v = R 2 . 

We discuss three possible limits of w e whose homogenized equation 
can be found. Namely, 

(11) u* = limsupi^, w* = liminf w e/ and u = limu £i , 

e^O 

where the subsequence u e . is such that E(w £l .) = \ j n \Vu £i \ 2 dx — > 
liminf^^o E(u E ) as i — > oo (See Figure 2). It turns out that w*, w*, 
and m don't follow simple homogenization whose boundary data is 
a simple average g. It means there are nontrivial homogenization 
processes for each different limits. For R ± 0, set 

~g (R, ei) = max g, g,(R, ^) = min g, 



and 

(12) ~g(R,M) = 



g\R,e x ) forM>supg, 
gXKei) forAKinfg, 
M for inf g < M < sup g. 



IfR = 0, let 



g(0^i) = ^(0 / ei) = g(0 / ^i) = g(0). 
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Proposition 7.1. For the particular choice v = e\, in equation (10), the 
limit functions u*, u*, and u will satisfy 



(13) 

where 

(14) 



A n = in Q, 

u(x) = M on x • e\ - -R lf 

u(x) -A onx- e\-R lf 



A = 



g\R 2 ,ei) ifu = u\ 
gA R 2,ei) ifu = u*, 
g{R 2 ,ei) ifu = u. 



Proof Select e z such that gij:) - max g and then u £i - u* since all u £i 
have the same boundary values. Hence it is clear that lim sup £ ^ u £ = 
u*. Similar argument can be applied to u* and u to have the conclu- 
sion. □ 



Proposition 7.2. In equation (10), when v + e\, u £ converges to u which 
is a solution to 



(15) 



A n = in Q, 

u(x) - M on x • v - -Ri, 

u ( x ) -g onx - v - R 2 . 



Proof Choose any point x such that x • v = R 2 and | = x (mod 1) for 
some x e Qi. Then, for v ei, 



(16) 



T g(-\da y -e 1 -vj g(^-)—-—dy 1 

J{(y-x).v=0}nQ r \ £ ' J { | yi _^|< r) \£/ei-V 

= f g\-r\ d y^^g> (as e ^ 0), 



which is independent of the choice of x. It implies g(x, v) = ^(x, v) = 
g(x). Let P(x - y) be a Poisson Kernel of the the domain Q. Then 

u £ (x) = M I P(x - y)doy + f P(x- y)g[—)da y . 

Jx-v=-Ri Jx-v=R 2 V £ / 

Letting e tend to zero and using first Lemma |3.3| and then Lemma 



3.2 we shall have 



lim I P(x-y)g(— )do y = [ P(x - y)gda y , 



'x-v=R 2 

which implies the conclusion 
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v = e x 



When v = e. 





Hill 




v * e, 



(a) figZa 



(b)fig2.b 



Figure 1. /zg-2.a s/zoz^s £/ze boundary x-v - R 2 when v - e\ 
or v ± e\. fig 2.b shows three possible limits (u*, u* and u) of 
u e . 

The next interesting question is to find the limit equation for the 
general converging sequence. In the following lemma, we will show 
there is a converging subsequence whose limit takes any value be- 
tween the supremum of g and its infimum oni-ei = R 2 . 

Proposition 7.3. For any A such that min g < A < max g, there is a 
sequence {u £i } converging uniformly to u such that its limit u satisfies 



Proof. There are Ej — > such that gyj) - A since g is 1-dimensional. 
Therefore all u £i have the same boundary values, implying u £i - 
u(x). □ 

7.2. Example 2. In this example we confine ourselves to R 2 . We con- 
sider the case when g is periodic only in X\ -direction and the domain 
Q is convex with two parallel flat parts of boundaries, orthogonal to 
v - e\. For exactness we consider the following stadium like domain 



(17) 



An = in Q, 
< u(x) - M on x • e { - -R lf 
u(x) - A on x - e\ - R 2 . 
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Let g(x) - g(xi) be 1-periodic, and independent of ^-direction. Now 
let u e be a solution of the following equation 

{Au e - in Q, 

Mx) = g(*) on da 

then u* and u* (see ( [TT] )) are sub- and super-solutions of ( |T8l ) respec- 
tively. In general, they are not solutions. 

In the next result we state that the homogenized boundary data 
may not be continuous even though g(x) is smooth. 

Proposition 7.4. There is a a smooth 1-periodic function g(x\) of one 
variable X\ and a subsequence of solutions {u £i } to equation ( [28] ), converging 
to u such that u is not continuous on dd. 

Proof. If X\ ^ ±R, then x - (xi,x 2 ) e <?Q satisfies IDD condition and 
then any limit will satisfy the boundary condition u(x) - g. Now 

we select £ z so that g(j:) - maxg ^ g. Then u £i has a converging 
subsequence to u which will be discontinuous at x - (R, 1). □ 
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